Electromagnetic modes in dielectric equilateral triangle resonators 
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Resonant electromagnetic modes are analyzed inside a dielectric cavity of equilateral triangular 
cross section and refractive index n, surrounded by a uniform medium of refractive index n'. The field 
confinement is determined only under the requirements needed to maintain total internal reflection 
of the internal electromagnetic fields, matched to exponentially decaying evanescent waves outside 
the cavity. Two-dimensional electromagnetics is considered, with no dependence on the coordi- 
nate perpendicular to the cross section; hence, independent transverse electric (TE) and transverse 
magnetic (TM) polarizations are described separately. A linear combination of six plane waves is 
sufficient within the cavity, whose wavevectors are related by 120° rotations and whose phases are 
related by Fresnel reflection coefficients. Generally, the mode spectrum becomes sparse and the 
minimum mode frequency increases rapidly as the index ratio N = n/n' approaches 2. For specified 
quantum numbers and N, the TM modes are lower in frequency than the TE modes. Assuming the 
evanescent boundary waves escape at the triangle vertices, TE modes generally are found to have 
greater confinement of the fields inside the cavity and much higher quality factors than TM modes. 

PACS numbers: 41.20.-q, 42.25.-p, 42.25.Gy, 42.60.-v, 42.60.Da 
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INTRODUCTION 



Confinement of electromagnetic fields only by to- 
tal internal reflection (TIR) offers a path for design 
of simple optical resonators or lasers with low volume 
and low threshold power. To this end, various two- 
dimensional (2D) geometries have been considered for 
their different influence on the mode properties, includ- 
ing disks, 01 cleavedB and etched|^ semiconductor tri- 
angles, squares,!^ ||,^ and zeohte ALPO4-5 hexagonal 
crystals.!^ ^ The problem of equilateral triangle res- 
onators (ETR) considered here is particularly interesting, 
because the symmetry present allows for a considerable 
simplification of the analysis. Good estimates for the 
ETR mode wavef unctions, frequencies, and lifetimes (or 
quality factors) are developed here, based on a set of six 
plane waves within the resonator, matched to each other 
by Fresnel factors, producing exterior evanescent waves. 

Ghang et al.^ analyzed the ETR wavefunctions un- 
der the assumption of Dirichlet boundary conditions 
(DBG), using an exactly known solution^, for trans- 
verse magnetic (TM) polarization only. In the presence 
of a dielectric-dielectric boundary, however, the fields do 
not go to zero on the boundary, instead, the field con- 
finement can be provided by TIR, assuming a refractive 
index in the cavity, n, greater than that of the surround- 
ings, n'. Huang et aZ.[^ jT^, 12| made a more complete 
analysis for an ETR with dielectric boundary conditions, 
approximately matching interior standing wave fields to 
exterior evanescent waves. Mode frequencies and qual- 
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ity factors Q were estimated[^ from 1000 to as high as 
20000, using the flnite-difference time-domain technique 
(FDTD) and Fade approximation p3[ , with the highest 
Q values associated with TM polarization. On the other 
hand, typical Q values from 20 - 150 were measured in 
photoluminescence (PL) experiments |^ on GalnAs-InP 
ETRs with edges from 5-20 /im. Recently I analyzed 
the ETR modes using the DBG approximation , and 
assuming the escape of evanescent boundary waves at the 
triangle vertices as the primary decay mechanism, esti- 
mated typical Q's from 10 - 500, with the highest values 
for TE polarization. Both of these theoretical approaches 
involve approximations, thus, it is important to consider 
an alternative description of the modes and clarify how 
the polarization influences the mode lifetimes. 

The DBG approximation applied previously^, |lj] 
clearly does not describe the fields correctly at a 
dielectric-dielectric boundary. The goal here is to use 
the general knowledge of the modes from the DBG anal- 
ysis, but do the correct matching of plane wave fields 
inside the cavity with evanescent fields on the outside. 
This matching is accomplished by employing the Fres- 
nel reflection coefficients correctly for all the plane wave 
components present in the cavity. It is assumed that all 
the electromagnetic field components do not depend on a 
z-coordinate perpendicular to the triangular cross section 
(i.e., longitudinal wavevector kz = 0). Then Maxwell's 
equations and associated dielectric boundary conditions 
lead to separated problems for TM and TE modes. Each 
polarization is controlled by one component of the elec- 
tromagnetic field that must be continuous across the 
dielectric-dielectric boundary. For the TM modes, the 
controlling wavefunction is the longitudinal electric field, 
tp = Ez; for TE modes, it is the longitudinal magnetic 
field, 2p = Hz. 
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The exact solution for an ETR with Dirichlet bound- 
ary conditions is a superposition of six plane waves of 
equal strengths but different phases and wavevectors. 
The waves undergo a — tt phase shift when reflecting from 
the faces, as required such that the incident and reflected 
parts cancel exactly at the boundary. Any one of these 
waves, when followed through a sequence of reflections 
due to TIR, returns to its original direction after six re- 
flections (see Fig. 1 of Ref. p^ . 

The situation is similar for the ETR with dielec- 
tric boundary conditions or Maxwell boundary conditions 
(MBC), where the correct field matching for Maxwell's 
equations must be applied. For either TE or TM po- 
larization, I show that a combination of six plane waves 
is still needed within the cavity, related to each other by 
the reflections from the three cavity faces. For a resonant 
mode whose fields are TIR-confined within the cavity, all 
the waves must impinge on the faces at incident angles 
greater than the critical angle Oc, where 



sin0, = l. 



n' 



(1) 



The generalization here compared to the DBC problem, 
is that when reflected from the faces, the waves undergo 
phase shifts determined by unit-modulus Fresnel coeffi- 
cients, as well as being rotated in propagation direction. 
The goal here is to determine the correct wavevectors and 
complex amplitudes of these six plane waves such that a 
fully self-consistent wavefunction is determined. 

When incident on the faces from inside the cavity, the 
waves produce evanescent waves in the exterior region 
just outside the cavity. It is difficult to give an exact 
description of the exterior fields, however, an approx- 
imate description is possible provided the six interior 
plane waves are incident on the cavity faces at angles 
sufficiently greater than 6c- In that case, the penetra- 
tion depth into the exterior region d — \k'j_\~^, is much 
less than the cavity size or edge a. Then, there are only 
strong evanescent fields very close to the cavity, and a 
description of these based on the transmission amplitude 
via Fresnel's equations is appropriate. The mode life- 
time and Q are estimated based on the assumption that 
the evanescent boundary waves radiate to the cavity sur- 
roundings when they reach the triangle vertices. This is 
a "strong damping" approximation, in the sense that the 
refiection or radiation of boundary wave energy back into 
the cavity is assumed to be insignificant. 

The presentation proceeds as follows. First, the MBC 
field matching of the six plane waves is described. Equa- 
tions for allowed quantum indexes are found, giving solu- 
tions for the mode wavevectors and wavefunctions. The 
wavefunction description is general enough that it de- 
scribes equally well the DBC, TM and TE problems, 
or any other boundary condition whose reflection phase 
shift is determined by incident angle. Lifetimes will be 
estimated using the boundary wave radiation, and com- 
pared to the simpler DBC theory and other analysis. 



II. EM FIELDS DESCRIPTION 
A. Six interior plane waves 

Here I describe the notation for the interior waves (ip), 
with assumed frequency to = c*k, where c* — cj ^feji is 
the light speed in the cavity, and the wavevector modu- 
lus fc is to be determined. The ETR has edge length a, 
and xy coordinates are used where the origin is set at its 
geometrical center, Fig. |l|. An initial wave @ is supposed 
to emerge from lower edge 6o at an angle a\ = ol rela- 
tive to the X-axis, and hence has wavevector components 
fci = (kx^ky) = (fccosa, fcsina). Only certain values of 
a will lead to a solution for an eigenmode of the cavity; 
determination of the possible values of a is an essential 
part of the solution presented here. 

This initial wave @ was formed by reflection of a wave 
incident on boundary bo at incident angle 9i^i = 90° — a. 
By consideration of the triangular geometry and using 
the law of reflection, it is seen that the incident angle 
when impinging on boundary bi will be 9i^2 — 60° — Oi^i, 
see Fig. As the wave subsequently reflects from the 
boundaries in ordered sequence bi, 62, ^o, etc., the inci- 
dent angles on each succeeding boundary simply oscillate 
between only two values symmetrically above and below 
60°. Thus, without loss of generality, one can assume 
a > 60° for the remaining analysis, and take k^ and ky 
of wave @ as both positive. The sequence of waves gen- 
erated from wave @, by sequential reflections from 61, 
^2, bo, bi and so on, are labeled as waves (2), (3), 0, @ 
and (6). Finally, wave (5) emerges from 62 and impinges 
on bo to regenerate an amplitude of wave @. To produce 
a consistent solution, this regeneration of wave Q must 
be in phase with the original wave @. 

Considerations of the law of reflection together with 
the equilateral geometry leads to the basic properties of 
the six waves, as summarized in Table The xy compo- 
nents of each wavevector are defined from the ai in the 
usual way, 



ki = A: (cos a/, sin a;). 



(2) 



It is seen that pairs of wavevectors, [A:i,fc6]j and [fc2,fc3] 
and [A:4,fc5], are related to each other by changing a — > 
—a. Each pair is related to the others by 120° counter- 
clockwise rotations around the z-axis, denoted by oper- 
ator R, which can be represented by the square matrix. 



R 



-1/2 -V3/2 
%/3/2 -1/2 



(3) 



The pair [fci,fc6] can be considered the original waves, 
from which the others are obtained. Wavevector fcg is 
simply obtained from fci by reflection across the x-axis. 

Due to the triangular symmetry, and the unit-modulus 
refiection coefficient under TIR, all six waves must have 
equal magnitudes, but different phases. Then the net 
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TABLE I: Definitions of the parameters of the plane waves 
labeled by wavevectors ki, within the triangular cavity, a; is 
the angle that each ki makes to the a;-axis. R is the operator 
for rotation through +120° around the £-axis. 



® 


® 


® 


@ 


© ® 


Q; = Of 


~a + 240° 


a + 240° 


~Q + 120° 


a + 120° -a 


ki — ki 


-ka 


R^ -ki 


R ■ fee 


R ■ ki ki {—a) 



wavefunction within the cavity is written as a sum over 
the waves 

6 

V'^J^Aie*"'-^^ (4) 
1=1 

where the Ai are all of unit modulus, and ip ~ for TE 
modes or tp = for TM modes. In both polarizations, 
■0 must be continuous across the cavity boundary. The 
TE and TM wavefunctions will not be the same, however, 
due to the different phase shifts implied by the Fresnel 
factors. 

Now consider what happens to these waves in terms of 
their interaction with the lower boundary, bo, which lies 
parallel to a; at y = — a/(2\/3). By appropriate symmetry 
transformations, the effects on the other two boundaries 
can be inferred. Three of the waves must be receding 
from bo and three must be approaching bg. Using the fact 
that 60° < a < 90°, for example, writing a = 60° + d, 
where d < 30°, the respective angles of the ki to the x- 
axis are 60° - z9, 180° - z?, -60° + z?, 60° - i9, 180° + d, 
and —60° — -d. Those that fall in the range from 0° to 
— 180° (including shifts by multiples of 360°) are incident 
on 6o. Thus, the waves @, (S) and © are incident on 
bo- Their incident angles are obtained from the formula, 
9i = 90° + ai, as —90° is the outward normal direction 
to bo- Then waves (5), (5) and @ have incident angles on 
bo, 

e,^3 = a - 30°, 6l,,5 = a - 150°, 0i^e = 90° - a. (5) 

One sees that 9i^e always has the smallest magnitude. 6i^5 
is less than zero because wave (5) has a negative wavevec- 
tor component along x- 

Extending such arguments to the other faces, the ray 
diagram of Fig. |l| summarizes the results, showing the 
incident angles of all the waves on each of the faces. It is 
important to notice that the only incident angles in the 
problem are the three of Eq. |^. For example, wave (1) is 
incident on 61, but at incident angle ^^^3. The symmetries 
represented in Fig. ^ will be exploited to determine the 
correct wave amplitudes. 

B. Field matching at the dielectric interface — 
Maxwell boundary conditions 

Now consider the reflections of the waves (5), (5), and © 
from bo, starting with wave ©. When wave © impinges 




L cx-60 bo 

FIG. 1: xy coordinates, with origin at triangle center, and 
boundaries bo, bi and 62. Rays represent the six planes waves 
propagating within the cavity, sketched for the value, a = 70° . 
The plane wave incident angles here determine the reflection 
amplitudes and phases used in Table 0. 



TABLE II: Relations between the incident and refiected 
wave amplitudes on the lower boundary (&o), the upper right 
boundary (bi) and the upper left boundary (62). The net re- 
flection phase shifts are Ai — Si — kiy-^, where 5i — S{9i^i) is 
the Fresnel reflection phase shift. See Fig. |^ for the geomet- 
rical reasoning behind this table. 
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incident 
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bo 
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@, Ai = Ase^-^^^ 


bo 


®, As 


6i,5 = Q — 


150° 


@, A2 = Ase'^^ 


bo 


®, Ad 


e,,6 = 90° 
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®, Ai = Aee'^'o 


bi 


®, A^ 






@, A2 = Aie'"'^ 


bi 


®, A3 


Oi,5 




®, Ae = ^36'^= 


bi 


®, A4 


di.6 




®, As = A4e*'^« 


b2 


®, As 


di,3 




®, Ae = Ase'""' 


b2 


®, Ai 


ei,B 




@, Ai = Aie'-^^ 


b2 


®, A2 






@, A3 = A2e''^« 



on bo, its reflection regenerates wave @ and it also pro- 
duces an exterior evanescent wave, denoted ©, of ampli- 
tude Fq (measured on the boundary) and wavevector fcg, 
with ke',x = fce.x by Snell's Law. Ignoring the other inte- 
rior waves, we match incident wave © and reflected wave 
@ to the exterior evanescent wave, ©. Matching the net 
-0 at y = -a/(2%/3), we have terms all proportional to 
the common phase factor, gHfe^a^-wt)^ where = kcosa, 

Age^'^^^-^^^ + Aie''''^-'-^^ = Fe- (6) 
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to 




71/2- 



71/4- 



FIG. 2; Dependence of the TIR Fresnel reflection phase shifts 
on incident angle, for TM [Eq. (H)] and TE [Eq. @] po- 
larization at the indicated refractive index ratios N, assuming 
unit magnetic permeabilities. 



On the other hand, in the usual analysis of the reflection, 
the correct boundary conditions for Maxwell's equations 
will relate the incident (■(/'inc.) and reflected (i/'rcfl.) ampli- 
tudes, measured at the boundary, by a Fresnel reflection 
factor r[9i), 



V'rcfl. 
V'inc. 



(7) 



Here 5(0i) is the phase shift experienced by the wave 
upon reflection, for incident angle 9i, in conditions of 
TIR. The 6'i-dependence of the different Fresnel factors 
appearing for TE and TM polarization is discussed sub- 
sequently in Sec. Ill, see Fig. for typical dependence 



with dielectric boundary conditions. 

Expressions ^ and imply that the wave ampli- 
tudes Aq and Ai have a relation. 



Axe 



Afi 



(8) 



The phase factor involves the y-component of k^, which 
is the negative of the y-component of ki. The subscripts 
on r and S indicate evaluating at the incident angle for 
wave ®, di,e = 90° — a. This allows us to write the 
reflected and evanescent field amplitudes in terms of Aq, 



A^^A^e^'^e-'"'-^^, ^6 = ^(1- 



73. (9) 



An equivalent algebra applies for analyzing how wave 
@ incident on produces reflected wave (2), and an 
exterior evanescent wave ©, in terms of incident angle 
^j,5 — a — 150°, and associated phase shift, 6^^. The re- 
flected and exterior field amplitudes are 



A2 - A5e'^=e-*'=^''73, ^5=^5(1 



(10) 



Finally, wave (3) incident on &o produces reflected wave 
and evanescent wave @, with similar expressions. 



A4 



i&3 



-ik3v-% 



v^, ^3=^3(1 



(11) 

The same type of analysis can be applied to the waves 
incident on boundaries 61 and 62, using the incident an- 
gles seen in Fig. ^ One can also use the symmetries under 
120° rotations R, seeing that, for example, R ■ bo = 61, 
and R ■ kg = k^, R ■ ki = fcs, implying that the phase 
relationship between waves (e) and (1) on bo is the same 
as the relationship between and on 62. Similar ar- 
guments apply to the other pairs of waves. The results 
of this analysis are summarized in Table where the 
phase factors relating the wave amplitudes are denoted 

by 

Ai=Si-kiy^, ; = 3,5,6. (12) 
v3 

There are only three distinct incident angles in the prob- 
lem, corresponding to the three fundamental Ot of Eq. ||. 
These factors A3, A5 and Ag are now the only quantities 
needed to correctly match self-consistently the interior 
fields, according to the nine equations in the last column 
of Table 0. 



III. DETERMINATION OF THE RESONANT 
FREQUENCIES AND WAVEFUNCTIONS 

The nine equations in the Table ^ can be summarized 
as three basic ratios. 



A4 
^3 



A2 ^ 
Ai ^ A5 



J^3 



A^ A^ Ar ' 



Ao Ai Ai 
Due to the simple structure, one sees the relation 
AiAaA^ 



12^4^6 _ ^i3A3 



(13a) 



(13b) 



(13c) 



(14) 



AxA-iA^ 

This leads to two fundamental relations for these phases, 

gi3(A3-HA6) ^ -^^ gj3(A5-l-A6) ^ Y. (15) 

Using all of the equations (Csl) together gives the relation. 



gi(A3+A5+2Ai3) ^ Y 



(16) 



Although these last three relations are not linearly inde- 
pendent, all are required to describe the solution. The 



5 



first two imply introduction of some integers denoted as 
ns and ns, such that 



A, 



2-K 



2tt 



Eq. ( p^ ) leads to, on the other hand, 

A3 + As + 2A6 ^ 27rn6, 



■715. 



(17) 



(18) 



where ng must also be an integer. Comparing these equa- 
tions demonstrates the constraint, 



"3 + "-5 = Sng, 



(19) 



that is, the sum of 713 and must be a multiple of 3. 
However, not all possible choices of these integers will 
lead to allowed solutions. Once the allowed quantum 
numbers and 715 are determined, the amplitude ratios 
of the six plane waves will be determined. 

Eqs. ( p7| ) and ( p^ merely give some sums of the A/ 
phase factors, whereas, we actually need to determine 
each one separately, and more importantly, we need to 
find the wavevector magnitude k. This is accomplished 
by using their definitions ( |l^ ) together with the cor- 
responding y-components of the wavevectors, obtained 
from (|). The necessary components are 

k 

ksy = fcsin(Q;— 120°) — — — (sina + VScosa), (20a) 



k sm{a + 120°) = (sin a - cos a), (20b) 



kyG = kasina 



71 f'""^ 



{S3 



266 



(23b) 



Remembering that the Fresnel phase shifts 5{9i,i) depend 
ultimately on a, via equations (^, the Eqs. ( p3| ) are seen 
to be coupled transcendental equations for unknowns k 
and a, assuming 713 and 715 are given. They can be solved 
in various ways. A simple approach is to eliminate k, and 
then determine the allowed a as the roots of the following 
function, 



P{a) 



2tt 

Y(»^3+"5)- ('53 + <55 + 2<56) 



%/3 



2-K 



(713 - ris) - (^3 - 5^) 



cos a (24) 
sin a — 0. 



In this last expression ng was eliminated using the con- 
straint (pj|). In the general case, it is not possible to solve 
for a in closed form. On the other hand, some straight- 
forward analysis of this function, together with numerical 
evaluation indicates the region in which to look for the 
quantum numbers 713 and 715. 

Once a has been determined, the modulus of the 
mode's wavevector, fc, can be found from either (|2^a) 
or (p^), or by their combination, written in a form like 
that for the known solution for the DBC case [Eq. (pq)]. 



ka = — { 
3 I 



773 - ^^5 



-^5) 



776 



1 

2^ 



{5^ 



2Sa 



-,2-, 1/2 
)] } -(25) 



key — fcsin(— a) = — fcsino;. (20c) 
Then this leads to the A; combinations, 

A3 + Ag = {S3 + Sq) + ifca(\/3sina -I- cosa), (21a) 



A5 -t- Ag = (^5 + Sq) + ^fca(V3sina — cosa), (21b) 



A3 + As + 2Ag = (^3 + (55 + 2(5g) + V3 /casina. (21c) 

The last of these then determines ky = k\y = k sin a. 
Another combination involves only kx = kix — fccosa, 

2tt 

A3 - A5 = {S3 - S5) + kacosa = -^{^3 - n^). (22) 

Then the basic wavevector components of wave (T) are 
expressed as 

27r 

kxa = kacosa = -^{^3 - ns) - {S3 - S5), (23a) 



A. The resonant wavefunctions 

Assuming a has been found (see below for DBC, TE 
and TM cases), then k^, k ^, a nd k or lo — c*k are deter- 
mined. In addition, Eqs. ( p^ ) determine the phase fac- 
tors needed to express the complete wavefunction for any 
mode, in terms of a via Eqs. (|2|) and (||) . The quantum 
numbers and Fresnel phase shifts give the A; as 



A3 = ^ 



A5 = - 
3 



2tt 

— {2n3-n5) + {S3 + S5-Se) 



2-K 

— (2n5 - 713) + {S3 + S^ - Sq) 



Ar = 



1 



27r 



-("•3 + "•5) - {S3 + (Js - Sq) 



all of which depend on a net phase factor, 
ip{a) = 53 + - (5g. 



(26a) 



(26b) 



(26c) 



(27) 
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This allows evaluation of the wavefunction (Q), simplified 
by freely choosing the amplitude of wave @ as 



where Aq is the overall wavefunction amplitude. 



ip = Aoi e*""^ cos 



2tt 



1 



9 



6' 



+ 6*''=*"'^ cos 



27r, „ ^ 1 " 
fcsj/y ~ Y^"^ ~ ^^^^ ~ 6*^ 



(28) 



, (29) 



Using values ns > and evaluating the possible quan- 
tum numbers, the well-known solutions for DBC are 
recoveredpl O, d3, in terms of shifted quantum numbers, 



TO = (713 — 715), n = {ne + 2). 



(34) 



These must be both odd, or both even, with the restric- 
tion m < n. Generally, we use these quantum numbers 
to label the modes with dielectric boundary conditions, 
in place of n^^nc,. Then the DBC wavevector moduli are 
reproduced from (I2F 



(35) 



where the wavevectors of waves ©, (5), and @ are needed, 
e.g., from Eqs. ([23|) and using Table |, 



1 



V3 



h 



«3y 



2 



(30) 



Note that the prohibited solutions at a = 60° would 
correspond to the prohibited case, m ~ n. It is straight- 
forward to check that the (to < n) DBC wavefunctions 
vanish on all the boundaries. 
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1 

^5y — ~2'^x — 2^yj (31) 



Note that if a = 60°, then 0i 3 = Oi^ = 30°, making 
6^ = (^6, and hence, </? = (Js = — tt (grazing incidence for 
wave (5), TE or TM pol arizat ion). Putting a = 60° also 
3 (See 



gives ^3 



2nH 



[II B 



below). These values cause 
the wavefunction "0 to vanish; there are no modes with 
a = 60° for the DBC, TE or TM cases, although it is 
tempting to sketch such ray diagrams. 



B. Equilateral triangle cavity with Dirichlet 
boundary conditions 

It is interesting to check the validity of the six-wave 
analysis in the case of Dirichlet boundary conditions; this 
also helps to locate the allowed 713 and 715 for Maxwell 
boundary conditions. The net field on the boundary be- 
comes zero when the Fresnel phase shifts are all taken to 
he 5i = —n (this is the limiting phase shift for grazing 
incidence, in either TE or TM polarization). Then the 
wavevector components reduce to 



kx Oi 



27r 



(n3 - ns). 



kya 



2-K 

71 



("6 ■ 



(32) 



The original assumption, a > 60°, imposes the relation 
tana = ky/kx > \/3. At the limiting value a = 60°, 
however, the net wavefunction ?/; of Eq. ( p9| ) vanishes; 
a = 60° is not allowed. There results < 2n5 -I- 3 
(equality maps to a = 60°). Also it was implicitly as- 
sumed that kx > 0, hence, we require > n^. There- 
fore, a given choice of ^5 allows a limited range of 77.3, 



n5 < < 2n5 + 3. 



(33) 



C. Maxwell boundary conditions: TM or TE 
polarization 

To define Eq. ( p4| ) for p{a) = 0, we require the Fresnel 
phase shifts for MBC. For TM polarization, tp = E^, and 
the Fresnel reflection coefficient for electric field polarized 
perpendicular to the plane of incidence is needed. For 
incident angle 9i, the formula isllF 



V'inc. 




(TM). (36) 



Unprimed electric permittivity e and magnetic perme- 
ability fi correspond to the cavity medium, whereas the 
primed values are those outside the cavity; the refractive 
indexes result from n = ^/ejl. Angle 9' is the refraction 
angle, a complex quantity obtained from Snell's Law un- 
der TIR, 



17' = i\J (sin6'i/ sin^c)^ 



1, 



(37) 



where critical angle 9c is defined in Eq. (|i|). The phase 
shift can be expressed also via 



6 /i / cos-^ ( 



(TM). 



(38) 



In the usual case for many optical materials with fi sa 
/i' « 1, the phase shift changes slowly as 9i ranges from 
9c {S = 0) to 90° {5 — —tt). Typical examples for index 
ratio N — 2.5, 4.0, 10.0 are shown in Fig. ^. In numerical 
results presented here we assume /i = = 1 , as for most 
optical materials. 

For TE polarization, tfj = H^, and the Fresnel refiec- 
tion coefficient for electric field polarized within the plane 
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80 a 85 

c 



FIG. 3: Example p{a) functions [Eq. @] at N = 8, = /i', 
for TM polarization, giving the lowest 13 valid mode solu- 
tions, located by zero crossings (dotted line). The individual 
curves are labeled by (m, n), see Eq. (p4[), where m changes 
in unit increments. There are three different groupings of 
curves. The bottom group (7 curves) has n — m — 2, the mid- 
dle group has n — m = 4, and the top group has n — m — 6. 
The values at the left limit are p(60°) — {n ~ m)n. The right 
limiting point is Oc = 90° — 6c- 



of incidence is needed. The Fresnel formula isvAi 



•0rcfl. 
V'inc. 



h cos Ui 



^ cos 6i + 



(TE). (39) 



Equivalently, the reflection phase shift is given by 



(TE). 




(40) 



The important difference, in comparison with that for 
TM polarization, is the presence of the factor e/e' instead 
of This significantly enhances the initial rate at 

which 5 increases with 9i, as can be seen in Fig. ||. Both 
the TM and TE phases shifts reach -tt at 0i = 90°, but 
for TE the approach is much more rapid. This can also be 
taken to imply that a DBC approximation for the modes, 
such as that applied in Ref . is more reasonable for TE 
modes than for TM modes. 

















or, 


Cm — M 




_ N=2.2 \ 










2.5\ 












3.2\ 












4.0\ 












5.0^ ^\ 






1 , , 




8.0^ 





60 



65 



70 



a 



75 



80 



85 



FIG. 4: Example p{a) functions for TE polarization, with 
the quantum numbers fixed at their lowest nonzero values, 
for various indicated index ratios N. Only the curves which 
cross zero (dotted line) give valid mode solutions (N > 4). 
The different termination points of the curves occur at Qc = 
90° - Be. 
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FIG. 5: The lowest mode wavevectors for TM polarization 
at index ratio N = 3.2, plotted versus quantum index m = 
^13 — ns. for different values of the other index n = n& -\- 
2 (Eq. indicated next to the curves, changing by unit 
increments. The solid circles indicate allowed modes; the solid 
lines connect those having equal values of n. The dotted lines 
locate the limits of stable TIR as expected from using the 
DBC solution for the ETR, Eqs. (^) and (0), explained in 
the text. 



D. Mode solutions: determination of a and (113,715) 
pairs, Maxwell boundary conditions 

Now consider the solution of Eq. ( p^ ) for a under 
MBC. As shown earlier, there are no solutions with 
a = 60°, because such a choice causes the wavefunction 
( p9| ) to vanish. The other extreme, a = 90°, which re- 
quires ria = 77,5, is prohibited because wave © can never 
experience TIR at a vanishing incident angle. Then, for 



MBC, similar to Eq. (js^), the search for possible (723, 715) 
pairs must take place in the range 

715 < 713 < 2715 + 3. (41) 

It is straightforward to calculate p(oi) [Eq. (|2j)] numer- 
ically and determine the zero crossings, which only take 
place provided that 713 is fairly close to the upper limit 
of (PI). Choices of 713,71.5 were made as follows. Start- 
ing from some 715 > 0, calculate 713 = 2715 -|- 3, which 
gives a prohibited pair at a = 60°. Automatically the 
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sum (na + ns) is a multiple of 3. Then the first pair 
to check for a valid solution, is to take reduced by 
1 ("3 ^ "-3 ~ 1) and increased by 1 {n^ + 1), 

such that the sum is the same multiple of 3. The result- 
ing pair (71.3,7^5) will likely have a solution for a if 715 is 
adequately large. This initial pair is equivalently set by 



713 = 2^5, 



ns = 1,2,3.. 



(42) 



Other pairs to try are found by continuing the reduction 
of by 1 together with the incrementing of ris by 1. 
The solution for a must occur within the range 



60° < a < 90° - 



(43) 



because the incident angle of wave ©, — 90° — a, 
cannot surpass the critical angle. Therefore, there is a 
considerably larger range for possible solutions for a as 
the index ratio increases (smaller 9c)- Conversely, an 
index ratio only slightly above 2.0 makes a substantially 
limited search range for a, as a result, fairly large values 
of ns are required before a solution is found. 

Plots of p{a) are surprisingly close to linear, especially 
for larger index ratio N and TM polarization. Examples 
are shown in Fig. ^ for index ratio N = 8, with various 
(^3, ns) pairs [or in terms of (to, n)] and TM polarization. 
One should note that the various p{a) curves at different 
quantum numbers are closely related; the obvious visible 
change occurs in the typical slopes as and 71-5 are ad- 
justed, while the Si{a) do not change. The dependence 
of some p{a) curves on N is indicated in Fig. ^, at fixed 
(7x3,715) = (2,1) for TE polarization. The terminating 
points of these curves occur at Uc — 90° — 9c- Then, the 
possibility for a zero crossing is greatly enhanced as N 
increases. Similar plots for TM polarization show curves 
much closer to linear form. The strong downward curva- 
ture as a —!■ etc for TE polarization means that for given 
quantum numbers, a TE polarized mode can be excited 
at lower index ratio than the same TM mode. 



IV. CALCULATIONS OF MODE SPECTRA 
AND PROPERTIES 

Having found (713,775) pairs and associated a, we can 
look at the mode dependence on polarization and refrac- 
tive index ratio N. Henceforth, modes will be labeled by 
the quantum number pairs, (771, 7i). 



A. TM polarization 

For TM polarization, the frequencies of the lowest 
modes are displayed in Fig. || for N = 3.2, in terms of 
their dependences on the mode indexes (777,71), defined 
in Eq. (|3J). This value for index ratio was used in Ref. 

in numerical and experimental studies of ETR semi- 
conductor cavities, with a different theoretical analysis of 
the modes. In general, the trend is for ka to increase with 
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FIG. 6: The lowest mode wavevectors for TM polarization 
at index ratio N = 8.0, plotted as described in Fig. |^. Note 
also that the fundamental mode here has lower (m, n) and ka 
than that for N = 3.2 . 



increasing quantum numbers. Just as in the DEC anal- 
ysis, the allowed modes for MBC must have indexes to 
and 77. either both odd, or both even (parity constraint). 
Note, however, that the TM mode wavevectors for cho- 
sen (to, 71) are noticeably lower than the values expected 
from the simplified DEC theory, expression (^5|), see Figs. 
I and |l| below. This can be attributed to the fact that 
the reflection phase shifts [from (^8|)] due not tend to- 
wards — TT even at very large index ratio N. 

The positions of the ka upper limits (see dotted lines 
in Fig. 1^) can be estimated quite accurately by using the 
TIR cutoff for a mode from the DEC theory|14|. Using 
the DEC solutions, all six plane waves can maintain TIR 
at all the boundaries only when 



TO > 77 



N2 - 1 



(44) 



Combining this with Eq. (35) for ka leads to the DBC- 
TIR limiting curve, expected to hold reasonably well at 
higher index ratio. 



uja 27rN 
ka = — < TO. 

c* 3 



(45) 



Additionally, the right limiting points of each curve cor- 
respond to quantum index pairs with 77 = TO-f 2, or equiv- 
alently, 773 = 2775. Then also using the DEC wavevectors 
(^), together with constraint to < 77, valid for either 
DEC or MEC, leads to the ka lower limit, 



ka 



Loa Att 

> TO. 

c* 3 



(46) 



Indeed, all the results for Maxwell boundary conditions 
lie between these results, the dotted lines in Fig. ||, lend- 
ing support to general aspects of the DEC theory. 

As N is increased, the lower TIR limit does not change, 
while the DEC-TIR upper limit moves steeper, or to the 
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FIG. 7: Fundamental TM modes of oscillation at N = 3.2, 
with (m,n) — (3, 5), where the pixel intensity is proportional 
to |7?e{?/'}|^''^, which enhances the definition of the nodal 
curves. The black nodal curves separate alternating regions 
of positive and negative Re{^}. Two degenerate wavefunc- 
tions are displayed. In a) phase is 9o — 0, in b) the phase is 
do ~ i"/2, where the mode amplitude is Aq = e'^", see text. 



upper left, encompassing more of the possible modes from 
the DBC theory, down towards smaller m for a given n. 
An example of this is given in Fig. ^, showing the corre- 
sponding results for index ratio N = 8. A larger number 
of states appears for each n, although the mode wavevec- 
tors have changed slightly. Q Conversely, as N decreases 
below the value 2.0, the DBC-TIR upper limit passes the 
lower limit, leaving no modes that can be confined by 
TIR. Furthermore, as long as N > 2.0, the total num- 
ber of modes is not finite, since n can be adjusted to an 
adequately large value to reach the fundamental mode. 
The main effect of placing N very close to 2.0 will be to 
force the lowest frequency mode to a large value of ka 



FIG. 8: First excited TM modes of oscillation at N = 3.2, 
with {m,n) — (4,6), as described in Fig. ^. Two degenerate 
wavefunctions are displayed. In a) phase is 6o = 0, in b) the 
phase is — tv/2, where the mode amplitude is Aq = e**", 
see text. 



and associated large minimum value of n. 

The fundamental mode quantum indexes depend on N . 
At N = 3.2, the fundamental mode has {m,n) = (3,5), 
and a diagram of its interior wavefunction (ip — E^) is 
shown in Fig. ^ The intensity of the pixels in these im- 
ages has been set proportional to |Re{^/;}|^/^, rather than 
linear in Ke{ip}, in order to sharpen the appearance of 
the zero crossings. The resulting nodal curves separate 
neighboring positive/negative regions of the wavefunc- 
tions. The calculation of mode lifetimes and quality fac- 
tors indicated on the wavefunction diagrams is described 
later in Sec. IV C . 



Wavefunctions for the first excited TM state at N = 
3.2, with (m,n) = (4,6), are shown in Fig. ^. It is im- 
portant to note that for the TM modes, the fields have 
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FIG. 9: TM mode wavevectors as functions of the index ratio 
N, for modes indicated by quantum indexes (m,n). 
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FIG. 10: The lowest mode wavevectors for TE polarization 
at index ratio N = 3.2, plotted as described in Fig. ^ The 
fundamental mode here has lower (m, n) and ka than that 
for TM polarization. On the other hand, at fixed {m,n), the 
mode wavevectors here are higher than those for TM polar- 
ization. 



substantial nonzero amplitudes even at the cavity edges. 
Furthermore, all these TIR-confined modes are doubly 
degenerate, since one can choose the two values ±kx and 
form degenerate pairs of states. Alternatively, different 
degenerate wavefunctions can be obtained by putting the 
overall amplitude Aq — e'^", choosing some arbitrary 
phase 00, and using only the real part of ([29|). The de- 
generate pairs presented here were formed using 9q = 
and 9o = tt/2. 

A rather weak dependence of some of the TM,„,„ mode 
frequencies on index ratio is shown in Fig. ^. The cut-off 
index ratios appear clearly as the left termination points 
of each curve. These cutoffs are similar in magnitude to 
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k.i = 12,5628 
zc'/a = 0.125 
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FIG. 11: Fundamental TE modes of oscillation at N — 3.2, 
with {m,n) — (2,4), as described in Fig. ^. Two degenerate 
wavefunctions are displayed. In a) phase is 9o = 0, in b) the 
phase is 9o ~ 7r/2, where the mode amplitude is Aq = e**", 
see text. 



that from the DBC theory, rewriting Eq. 



N > Ne= \I3^ + 1. 



(47) 



(See Fig. |l| for the DBC cutoffs.) However, the TM 
mode wavevectors are considerably lower than the pre- 
diction of the DBC theory, Eq. (|3q), because the TM 
reflection phase shifts 6i^i are never near — tt. 



B. TE polarization 

Results for the TE mode spectrum at N = 3.2 are 
shown in Fig. m], in the same manner as displayed for TM 
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FIG. 12: The lowest mode wavevectors for TE polarization 
at index ratio N = 2.5, plotted as described in Fig. H. 



polarization. Again, the allowed mode indexes satisfy the 
parity constraint, and the ka values fall within the limits 
of the DBC-TIR theory. On the other hand, there are two 
primary differences compared to TM polarization. First, 
the fundamental mode has a lower (m, n) pair, and in 
addition, a slightly lower ka than for the fundamental 
mode with TM polarization. Second, at a given (m,n), 
we see that the wavevectors for TE polarization always 
are noticeably higher than for TM polarization. 

At N = 3.2, the TE fundamental is (m, n) = (2,4), 
Fig. however, the Q for this mode is extremely low, 
so that it cannot be considered stable. 

In Fig. |l2| we illustrate the effect of decreasing N down 
to the value 2.5, closer to the extreme limit 2.0 . Two 
low modes at m = 4 and m = 5 have been squeezed 
out by the lower DBC-TIR upper limit, as well as the 
entire spectrum becoming narrower. Of course, a similar 
squeezing effect takes places for the TM spectrum. 

The dependence of some TE mode wavevectors on in- 
dex ratio is shown in Fig. |l^. Compared to the TM 
modes, the TE modes show a stronger variation with 
N, especially just above the cutoff ratio. The plot also 
shows the DEC mode wavevectors as dotted lines, termi- 
nating at the cutoffs predicted by the DBC-TIR theory, 
Eq. (^). Once N reaches adequately large values, the 
wavevectors from the calculations for dielectric bound- 
ary conditions asymptotically approach the DBC values. 
This can be attributed to the extra factor of e/e' in the 
TE phase shift formula (|40|), which easily causes all the 
reflection phase shifts to rapidly approach — tt, the value 
under Dirichlet boundary conditions. 

The fields diminish more rapidly near the cavity edges 
for TE polarization than for TM polarization. At N = 
8.0, the fundamental modes have {m,n) — (1,3), whose 
corresponding wavefunction is the simplest possible, hav- 
ing a single nodal line across the cavity, Fig. |l^. The 
avoidance of the fields close to the cavity boundaries of 
the TE mode is apparent, whereas, the TM fields tend 
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FIG. 13: TE mode wavevectors (solid curves) as functions of 
the index ratio N, for modes indicated by quantum indexes 
(m,n). The dotted lines show the DBC wavevectors, termi- 
nating at the cutoffs given by Eq. (^7[). 



towards maximum values near the boundaries. 

Taken together, these results emphatically confirm the 
idea presented in Ref. 
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that the simple DBC theory is 
much more appropriate for TE polarization than for TM 
polarization. 



C. Estimates of mode lifetimes due to boundary 
wave emission 

The (doubly-degenerate) solutions found here are ap- 
proximate. The assumption of evanescent traveling 
waves along each boundary, rotating around the trian- 
gle in the exterior region, is not exact, because the edges 
are of finite length. In reality, any of the evanescent 
waves can be expected to scatter from the triangle ver- 
tices, leading to radiation away from the cavity, and re- 
flection of that evanescent wave backwards from the ver- 
tices. There should be linear combinations of evanescent 
traveling waves moving in both directions along the edges 
(standing waves). This in turn would lead to a mixing 
of the doubly degenerate modes, i.e., the degeneracy will 
be split due to the scattering experienced by the fields at 
the triangular vertices. 

In the case of modes whose exterior wavelength is small 
compared to the cavity edge, these effects may be small, 
and ignoring the resulting degeneracy splitting, we can 
try to estimate the power loss, and hence, the mode life- 
time. To accomplish this, following Wiersig[p^, it is as- 
sumed that all the power in the evanescent boundary 
waves is radiated when reaching the vertices of the trian- 
gle, without reflecting backwards from the vertices. The 
lifetime r and Q-factor are calculated by 



U/P, 



Q 



27r/r, 



(48) 



where U is the total energy in the cavity fields, and P is 
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FIG. 15; Mode quality factors for a) TE and b) TM modes 
calculated from boundary wave emission, see (js^) and dss] ) , as 
functions of the index ratio N, with modes labeled by (m,n). 



FIG. 14: Fundamental a) TE and b) TM modes of oscillation 
at N = 8.0, with (m,n) = (1,3), as described in Fig. ^ both 
for phase 61 = 0, see text. (The choice d = 7r/2 would instead 
produce a vertical nodal line.) 



the total power radiated, as found from leakage of the 
evanescent boundary waves at the vertices. This ap- 
proach was used also to estimate the mode lifetimes using 
the DBC solutions jljl; there it was found that generally 
speaking, the TE modes have longer lifetimes than the 
corresponding TM modes, provided the index ratio N is 
large. It is important to check the relative mode lifetimes 
using the more correct dielectric boundary conditions. As 
the present solution has demonstrated the avoidance of 
TE fields near the cavity boundaries, one might expect 
the TE lifetimes to be longer. 

Following the calculations in Refs. |lj, |l^, the total 
energy of the fields in a cavity of height h can be written 



U = [ hdxdy = [ hdxdy ^^ffl; (49) 

J on J 8tt 

the two forms convenient for TM modes {\E\ = =11^) 
and TE modes {\H\ = = respectively. Using 
the wavefunction of Eq. (p9|), it is not possible to find 
a simple formula for this integral in the general case, 
thus it was evaluated by numerical integration within 
the triangle. In the DBC limit, with (p — —it, one finds 
/ ^l^liZ-p = 3\/3|Ao|V8; generally also for MBC the 
integral is of this order. 

The total boundary wave power is taken as the sum 
of the powers from boundary waves on each edge. At 
each edge there are three different incident waves, each of 
which generates an evanescent wave. For example, on bo, 
waves (5), (s) and @ separately produce the evanescent 
waves @, © and ©. We sum the powers in each of 
these evanescent waves, then multiply by three due to 
triangular symmetry, to give the total power radiated. 

The boundary wave power of an individual evanescent 
wave is calculated using the Foynting flux along the ex- 
terior side, S' = ^Re{i?' x H'*}, expressed alternately 
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as 



S' 



\E 



87re' 



\H 



(50) 



where x points parallel to the boundary, and the incident 
angle from within the cavity is 9i. The evanescent wave 
has an exponentially decaying behavior into the exterior 
medium. If Ai is the amplitude on the interior side, with 
boundary at y = 0, and the y-coordinate points from the 
boundary into medium n', Eqs. (0) and (p^ give 



iSiOi 



k'i'y\2 



(51) 



Integrating the total flux contained from y ^ Q to y ~ oo 
gives the power of this wave along the boundary, Px = 
h dy S ' ■ X, with only a tiny formal difference for the 
two polarizations. 



IA;Pcos2fi 



^1- (sine,/ sine,)' 



(TM), 
(TE). 



(52) 



In practice, however, the phase shifts for TE polarization 
are typically much closer to — tt than for TM polarization, 
which causes the TE powers to be smaller. The incident 
squared wave amplitude is l^ip = ;||Aop, as all the six 
wave components of ip are of equal strength, Eq. (^). 

The total emitted power from all edges is three times 
the sum of the powers on edge bo. 



P = 3(P,^3 + Pcc,5 + Pxfi). 



(53) 



Taking the net energy/power ratio and simplifying leads 
to the dimensionless lifetime expressions. 



TC 



1^0 



3,5,6 



-yi-(sin 0J sin Oi)^ 



t (54) 
^ (TE). 



The factor of e'/e < 1 for TE polarization tends to reduce 
the lifetime, however, the phase shifts in the denominator 
have an even larger effect, such that typically, the TE„i,„ 
lifetime is found to be longer than the TMm,n lifetime. 
The summation of power terms in the denominator is 
usually dominated by that of wave ©, which has the 
smallest incident angle. The expression is scaled by a/c*, 
the time for the signal to cross the cavity. A well-defined 
mode should have rc*/a ^ 1, however, usually it is more 
typical to look at the related mode quality factor, Q, 
defined by 



Q — LOT — c*kT — {ka){Tc* /a), 



(55) 



whihc is the dimensionless mode wavevector times the 
dimensionless lifetime. 

Numerical results for Q versus index ratio are displayed 
for the lowest modes in Fig. for TE and TM polar- 
izations, assuming ^ = fi' . One finds that Q's for TE 
modes are considerably larger than for TM modes, espe- 
cially far enough above the cutoff N for a given mode. 
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FIG. 16: Comparison of TEs.y wavefunctions at a) N = 2.5, 
only slightly above cutoff, and very weakly bound with small 
Q and b) N = 8.0, substantially above the cutoff, with much 
higher Q, and field amplitude concentrated away from the 
boundaries. 



Furthermore, the TM quality factors tend to saturate at 
large N, while the TE quality factors tend to increase 
proportional to = e/e'. This latter effect can be seen 
due to the asymptotics for the TE reflection phase shift, 
based on the identity. 



cos^ 6i 



sin^ 0c + (eVe)2cos2 Oi 



(56) 



The net result is that Q and r for TE modes are propor- 
tional to e/e' in the limit of large index ratio, a factor 
not present for TM polarization. 

It is interesting to observe the changes in the wave- 
functions with increasing N or effectively, increasing Q. 
In Fig. |l^, TE5J mode wavefunctions are plotted for 
N = 2.5, just slightly above the cutoff index ratio, and 
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FIG. 17: Comparison of some mode a) lifetimes and b) qual- 
ity factors, as functions of the index ratio N, from the present 
theory, using dielectric boundary conditions (solid curves), 
and from the simpler theory using Dirichlet boundary condi- 
tions (dashed curves). The modes are labeled by (m, n). 



for N = 8, well above the cutoff. At lower N (and Q), a 
central lobe of V' is connected to one vertex of the bound- 
ary, and another lobe is nearly coupled to the opposite 
edge. At higher N (and much higher Q), these interior 
lobes have become completely detached from the bound- 
ary, and the fields appear more concentrated within the 
interior, farther from the boundaries. 

A comparison can be made between the results using 
dielectric boundary conditions and those from the sim- 
plified DBC theory 0. Fig. |l^ shows Q and t results for 
the modes (1,3) and (2,4), as derived from the two ap- 
proaches, using the total power in all boundary waves for 
this comparison. Indeed, there are only minor differences 
between the lifetimes, mostly near the cutoff index ratios, 
due to the fact that the DBC theory overestimates the 
cutoffs. For the quality factors, the discrepancy between 
the approaches is much greater for the TM modes, pri- 
marily because the DBC theory has overestimated their 
frequencies (or wavevectors) . Nevertheless, the ability of 
the simplified DBC theory to predict the general trends 
with N is remarkable. 



TABLE III: Properties of some of the lower TM mode for 
N = 3.2, labeled by {m,n) or by underlined indexes (m, Z), 
where compared to results of Ref. in parenthesis. 



mode 


a 


ka 


Q 


TM3,5 


70.4167° 


14.3406 


32.18 


X iVl A (i T A\ 


68.6412° 


18.5733 (18.8) 


67.67 


1 IV 1 K 1 f n 1 ^ 


67.3705° 


22.7686 (22.9) 


107.9 


TMg^ (0,18) 


66.4225° 


26.9545 (27.0) 


154.5 


TM7^g (0,20) 


65.6894° 


31.1374 (31.2) 


207.9 


-L iVic in (Ti 991 


65.1060° 


35.3196 (35.3) 


268.4 


1 JVif; in /"I 911 


70.6463° 


33.5168 (34.1) 


145.8 




64.6307° 


39.5019 (39.5) 


335.9 


TM? 11 n 9^1 


69.6339° 


37.7112 (38.3) 


227.4 


TMio,i2 


64.2363° 


43.6845 


410.4 


TM8,12 


68.7916° 


41.8887 


307.3 


1 VI 1 c 1 '7 CM o^;^ 


62.9701° 


64.6045 (64.5) 


889.9 (15000) 


TMi3^17 (1,35) 


66.1053° 


62.7505 (62.9) 


784.1 (4380) 


TMii,i7 


69.3950° 


61.0455 


593.5 


TMig^lg (0,38) 


62.8025° 


68.7897 (68.7) 


1007 (22000) 


TMi4^18 (1,37) 


65.7527° 


66.9250 (67.0) 


899.0 (7000) 


TMi2,18 (2,36) 


68.8431° 


65.2085 (65.5) 


717.7 (1230) 


TMl7,19 


62.6528° 


72.9751 


1132 


TMl5,19 


65.4386° 


71.1006 


1021 


TMi3^19 (2,38) 


68.3516° 


69.3713 (69.5) 


845.1 (1960) 


TABLE IV: Properties of 


some of the lower TE modes for 


N = 3.2, labeled by (m,n) 


or by underlined indexes {rn,l), 


where compared to results of Ref. |lC|, in parenthesis. 


mode 


a 


ka 


Q 


TE2,4 


71.5207° 


12.5628 


1.576 


TE3,5 


69.9852° 


17.6896 


18.95 


TE4,6 


68.4450° 


22.1165 


60.13 


TE5_7 (0,16) 


67.2586° 


26.3873 (26.2) 


121.4 


TE6,8 


66.3498° 


30.6105 


200.1 


TE7,9 


65.6382° 


34.8149 


294.6 


TE5,9 


71.3190° 


32.0399 


15.63 


TEsao (0,22) 


65.0678° 


39.0107 (39.0) 


404.1 


TEe.io 


70.3771° 


36.7552 


59.80 


TEi4_i6 (0,34) 


63.1471° 


64.1404 (64.2) 


1350 (6130) 


TEi2,i6 (1,33) 


66.4654° 


62.2279 (62.4) 


789 (4100) 


TEio,i6 


69.8843° 


60.2288 


221 


TEi5,i7 (0,36 


62.9597° 


68.3271 (68.4) 


1560 (10860) 


TEi3,i7 (1,35) 


66.0736° 


66.4147 (66.6) 


972 (15320) 


TEii,i7 


69.2967° 


64.4910 


346 
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D. Comparison with other ETR theory for 
dielectric boundary conditions 

ETR modes have previously been analyzed by Huang 
et aL(HGW)pO|] using different approximations, also in- 
volving matching of interior fields undergoing TIR to ex- 
terior evanescent fields. In particular, HGW used what 
was called a "perfectly confined approximation for the 
transverse wavefunction." This is a Dirichlet boundary 
condition on part of the full wavefunction, giving the y- 
component of wavevector © [Eq. (9) of Ref . ^ , 



k^a = —=[m 
V3 



1), 



0,1,2, 



(57) 



Indeed, this is the same as Eq. ( |32[ ) reviewed here for 
the DBG problem, with m + 1 equivalent to rt = ng + 2. 
Thus, it is noticeably different from the result for dielec- 
tric boundary conditions, Eq. (]23|b), which more fully ac- 
counts for all the reflection phase shifts. The factor m + 1 
begins at the value 1, whereas, in our results, the mini- 
mum value of equivalent quantum index n is 3. For the 
other component of wavevector ©, the result obtained 
[Eq. (21) of Ref. |l§ is similar to Eq. (||a). 



20. 



1 = 3,4,5,, 



(58) 



where I has the same parity as m, and 6 is an indi- 
vidual TIR Fresnel phase shift, having different forms 
for TM and TE modes. Index I appears equivalent to 
m = — n^, however, Eq. (p3|a) involves phase shifts of 
two different waves, instead of the individual phase shift 
i- 

For a chosen value of n, different choices of m give 
solutions with ka of similar magnitude, hke the modes 
connected by solid lines in Figs. ||, |6[ |o| and |l^. This 
"number of transverse modes" increases with both n and 
N, and is smaller than that found in Ref. |lO[ Because 
ka = 27rN/A depends on (m,n) in a nontrivial way, it 
is not possible to give a simple expression for this mode 
count. 

For N = 3.2, a summary is made of the lowest TM 
and TE modes in Tables [II and IV and compared to 



results from HGW. In spite of the obvious differences in 
these theoretical approaches, they both lead to very sim- 
ilar predictions for the mode wavevectors or frequencies, 
agreeing to within about one percent. The prediction of 
the Q-values are considerably different; HGW used the 
finite-difference time-domain (FDTD) technique fl^ com- 
bined with Fade approximates to estimate Q. In partic- 
ular, the FDTD technique predicts that TE polarization 
results in much smaller Q than TM polarizationpH| , ex- 
actly opposite to our results (Fig. |l^). The short lifetime 
for TE modes was explained by the zero in the reflectivity 
at the Brewster angle Ob, which is only slightly less than 
the critical angle. However, in a confined TIR wavefunc- 
tion, all six plane wave components must impinge on the 
faces at angles greater than 9c and hence greater than 9b, 



TABLE V: Mode frequencies and free space wavelengths in 
the range around 1.3 to 1.6 ^m, for cavities with N — 3.2 and 
edge lengths a. The modes are labeled by m, n. Wavelengths 
marked with an asterisk fall about 2% below peaks in PL 
data, Ref. ^ 



a (//m) 


mode 


ka 


A (/im) 




2 


TMs 5 


14.341 


1.5675 


191.25 


2 


TEi 


12.563 


1.7893 


167.54 


5 




33.517 


1.6767 


178.79 


5 


TMo in 




1.5911* 


188.41 


K 
o 


_L iVi7^11 




1.4902* 


901 1 7 


5 


TMg.ii 


39.502 


1.4227* 


210.72 


5 


TM8,12 


41.889 


1.3416* 


223.46 


5 


TMio,i2 


43.685 


1.2865* 


233.03 


5 


TE5,9 


32.040 


1.7540 


170.92 


5 


TE7,g 


34.815 


1.6142 


185.72 


5 


TEe.io 


36.755 


1.5290 


196.07 


5 


TEs.io 


39.011 


1.4406 


208.10 


5 


TE7,ii 


41.135 


1.3662 


219.43 


5 


TEg.ii 


43.202 


1.3008 


230.46 



so it is hard to understand how the Brewster minimum 
comes into play, unless diffractive effects at the triangle 
vertices are strong, causing violation of the six-wave as- 
sumption. 



E. Mode generation near 1.55 /^m free space 
wavelength 

There has been recent interest, both experimental 
and theoretical, in using semiconductor ETRs operat- 
ing around 1.55 /i m wavelength. |l^, |ll|, |l^ Here we 
assume a cavity medium with n — 3.2, surrounded by 
vacuum, and summarize the mode spectra obtained by 
the present theory, for some typical cavity edge lengths, 
a — 2/im and a = 5/im. Of course, all mode frequencies 
simply scale as 1/a, thereby giving the primary tuning 
control parameter. Modes whose wavelength outside the 
cavity ranging from approximately 1.30/im to 1.60/im are 
considered; this corresponds to frequency range 187 THz 
< / < 231 THz. 

Results for all the modes found from 1.3 to 1.6 ^m are 
shown in Table ^ for TM and TE polarization. Obvi- 
ously very few modes occur in any such narrow range if 
the cavity is small (a = 2fim), which could allow for fine 
frequency tuning. Gonversely, many modes are present 
in larger cavities (such as a = 10/im) and single mode 
operation is difficult. A cavity with a = 5/zm has a 
moderate number of modes; the TM mode wavelengths 
are similar to those found in photoluminescence (PL) 
experiments!^. In fact, five of TM mode wavelengths 
calculated here are about 2% lower than five peaks seen 
in Fig. 2 of Ref. p|. This deviation might be attributed 
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primarily to an uncertainty in the cavity size, and sec- 
ondly, caused by a weak variation of dielectric constant 
with wavelength.]^ Discounting these factors, the agree- 
ment for the 5/im cavity is reasonable. Comparison with 
the experimental spectrum for a — lO^m is more dif- 
ficult, although it appears that some of the TMm_„i+2 
modes from the present theory (not shown) do appear 
in the PL data, again allowing for uncertainty in a and 
dispersion. 



V. CONCLUSIONS 

The phase relationships between the six plane waves 
within an ETR have been determined so they match cor- 
rectly to each other and to exterior evanescent waves, 
according to Fresnel reflection coefficients for a dielec- 
tric on dielectric boundary at index ratio N = n/n'. The 
theoretical wavefunction description is actually very gen- 
eral; it applies to any choice of reflection phase shifts, 
including those for DBC, and TE or TM polarizations. 
The main approximation here has been that the evanes- 
cent fields do not perturb the interior fields; assuming 
the evanescent fields radiate at the triangle vertices (a 
strong damping approximation), the mode lifetimes and 
quality factors have been estimated. Modes with high Q 
should be very well described by the wavefunction (29). 
The mode wavevectors are consistent with previous anal- 
yses [|l^, |ll|, 0, although the wavefunction description 
and quality factors differ from the FDTD technique ||Tl|] 
results. Some mode results are consistent with PL ex- 



perimental data on 5/im semiconductor ETRsQ. 

For N adequately above the cutoff for a mode (m, n), 
the TE mode wavevectors are very close to the predic- 
tions from the simplified DBC theory, whereas, the TM 
mode wavevectors are consistently below the DBC re- 
sults. This is because large N results in Fresnel reflec- 
tion phase shifts very near — tt, the value for the DBC 
theory, only for TE polarization. This causes the TE 
wavefunctions to avoid the cavity edges; on the other 
hand, the TM wavefunctions have significant amplitude 
at the edges. For both polarizations, dielectric boundary 
conditions give lower cutoff index ratios than from the 
DBC theory. As N increases starting from 2.0, the the- 
ory demonstrates how the spectrum of available modes 
expands (Fig. ||, etc.), but always stays within the limits 
predicted from the DBC theory. 

An extra factor of e/e' appears in the lifetime for TE 
modes, not present for TM modes. This theory then 
predicts considerably larger lifetimes and Q's for TE po- 
larization, whereas the FDTD approach |ll| led to larger 
Qs for TM modes. At large index ratio, the lifetimes 
found here approach the values found in the simpler DBC 
theory iQ, for both polarizations. 
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